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ABSTRACT 

We study the effects of stellar collisions, particularly on feeding massive black holes 
(BHs) and color gradients, in realistic galactic centers. We find that the mass released 
by stellar collisions is not sufficient to account for the present BH mass in galactic 
centers, especially in bright galaxies. This study, together with the study by Magorrian 
& Tremaine (1999) on tidal disruption of stars by massive BHs, implies that the 
material for BH growth (especially in galaxies brighter than ~ 10 9 Lq) can only come 
from other sources, for example, the mass released by stellar evolution in the initial 
~ 1 Gyr of the galaxy's lifetime, or the gas that sinks to the galactic center in a galaxy 
merger. We also analyze how the color of a stellar system is affected by collisions of 
stars. We find that collisions between main-sequence stars cannot cause observable 
color gradients in the visible bands at projected radius R > 0.1" in M31, M32 and 
other nearby galactic centers. This result is consistent with the lack of an observable 
color gradient in M32 at R > 0.1". At even smaller radii, the color differences caused 
by collisions between main-sequence stars are at most 0.08 mag at R — 0.02". The 
averaged blueing due to stellar collisions in the region R < 0.1" of M32 should not 
be larger than 0.06 mag in color index U — V and 0.02 mag in V — I. The observed 
blueing in the center of the galaxy M31 (in a 0.14" x 0.14" box) must be caused by 
some mechanism other than collisions between main-sequence stars. 

Key words: black hole physics - galaxies: individual: M31, M32 - galaxies: kine- 
matics and dynamics - galaxies: photometry - galaxies: nuclei - stars: blue stragglers 



1 INTRODUCTION 

The centers of galaxies are extreme astrophysical environments. They house massive black holes (BHs) (e.g., Magorrian et al. 
1998) and also densely packed stars (10 2 -10 4 Mq pc -3 at galactic radius ~ 10 pc, e.g., Faber et al. 1997). High stellar densities 
lead to frequent stellar collisions. The outcome of the collision between two stars depends on their types (e.g., evolutionary 
phase, mass and radius) and their kinematic parameters (e.g., relative velocity and impact parameter) (e.g., Freitag & Benz 
2002). The two stars may both survive their collision, or the collision may destroy one or both of them, or lead to their 
coalescence. Gas released from the colliding stars may be accreted onto the central BH, and this process has been proposed 
as one of the contributions to growth of central massive BHs (e.g., Frank 1978). Coalesced stars may appear as new stars 
with different luminosity properties or different colors from their parents; thus the luminosity and color of galactic centers 
are possibly affected by stellar collisions. The purpose of this paper is to study the effects of stellar collisions, particularly on 
feeding BHs and color gradients, in realistic galactic centers. 

Studies of central BHs in nearby galaxies have revealed a tight correlation between central BH mass and galactic velocity 
dispersion (e.g., Tremaine et al. 2002 and references therein), which strongly suggests a close link between the formation 
and evolution of galaxies and their central BHs. In an isolated stellar system, central BHs may accrete stellar mass from the 
following three sources: (i) gas released by tidal disruption of stars and stars swallowed whole by central BHs (mostly for 
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stars on elongated radial orbits which may come close to central BHs) , (ii) gas released by stellar collisions, (iii) and gas lost 
by stellar winds. Whether the above three sources are sufficient to account for the present BH mass and which is the main 
contributor to the BH growth depend on the stellar density distribution and velocity dispersion (e.g., Duncan & Shapiro 
1983; Murphy, Cohn & Durisen 1991; Freitag & Benz 2002). The first source in this list, tidal disruption of stars, has been 
studied by Magorrian & Tremaine (1999). They show that the total mass of stars that are tidally disrupted or swallowed 
whole over the lifetime of typical nearby galaxies is of the order of 10 6 M Q , approximately independent of galaxy luminosity. 
Thus disrupted stars may contribute significantly to the present BH mass only in faint galaxies (< 10 9 L Q ). One part of this 
paper (§ 3) studies whether the second source (i.e., gas released by stellar collisions) is sufficient to explain current BH masses 
in realistic galactic centers. 

The colors of galactic centers reflect the constituents of their stellar population, and radial color gradients may provide 
much information on the formation and evolution of the stellar populations or galaxies. In § 4, we will analyze how the color 
of a stellar system is affected by collisions of main-sequence stars. We will focus on studying two galaxies in the Local Group: 
M31 and M32, which are the nearest giant spiral and dwarf elliptical galaxy (distance: ~ 800 kpc) and have high-resolution 
observations of their centers in several color bands (Lauer et al. 1998). Multicolor HST WFPC2 (Wide Field Planetary 
Camera 2) images show that the center of M31 (the region within ~ 1 pc of its central BH) appears bluer than the rest of 
nucleus and the surrounding bulge. Lauer et al. (1998) suggest that this blueing is caused by collisions between main-sequence 
stars since coalescence of main-sequence stars by collisions may form blue stragglers (which are brighter and bluer than their 
parent stars, e.g., Leonard 1989; Bailyn & Pinsonneault 1995). However, observations show that the center of M32 lacks 
color gradients at similar spatial resolution. The lack of color gradients in M32 is inconsistent with a rough estimate of the 
effects of collisions between main-sequence stars by Lauer et al. (1998), and this inconsistency is claimed to be an important 
puzzle. In § 4, we will study whether collisions of main-sequence stars can cause observable color gradients in the centers of 
M31, M32 and some other nearby galaxies. 

Finally, our conclusions are given in § 5. 



2 STELLAR COLLISION RATES AND TIMESCALES 

We first present a general description of stellar collision rates and timescales. 

In a stellar system, the distribution function (DF) /(x, v) is defined so that /(x, v) d 3 xd 3 v is the number of stars 
within a phase-space volume d 3 xd 3 v of (x, v). We define the stellar mass function of the system at time t as S(m, t), so that 
H(m, t) dm is the probability of finding a star with mass in the range m — > m + dm (i.e., J °° H(m, t) dm = 1); and we define 
the stellar mass and radius function at time t as £(m,a,t), so that £(m, a, t) dmda is the probability of finding a star with 
mass and radius in the range (m, a) — > (m + dm, a + da) [i.e., J °° £(m, a, t) da — H(m, t)]. We assume that the stellar mass 
and radius function is independent of the location of the stars in the phase space (x, v). 

The cross section for a physical collision (or contact encounter) of two stars with mass and radius (mi,ai) and (m 3 ,a 3 ) 
moving with the relative velocity at infinity v rci is given by: 

V- I \ (2 2 (, . V ij\ _ / l^p if Vrcl > Vij 



where 6 is the impact parameter, r p = a t + aj, and — y^2G(mi + mj)/r p . If (m;, a^) = (m 3 , aj), = y /, 2Gm i /a i = v CBC is 
the escape velocity from the surface of either star. Thus, a star (mj, aj) moving with velocity Vj through a background of stars 
with distribution function /(r, v), suffers collisions with stars with mass and radius in the range (wij, aj) — > (mi + dmi, ai+dat) 
at a rate F i; , (r, Vj) dm^dai given by (cf., eq. 8-116 in Binney & Tremaine 1987): 

Y l3 (Y,v 3 ,t)dmidai = i(mi,ai,t)dmidai J "d 3 Vi /(r, Vj)|vj - v 3 |S i3 (|v, - v 3 1). (2) 
If the distribution of field stars is isotropic in velocity space, we set v = |v| and we have 

f°° 2-KVi r* +Vj 

Y lJ (Y,w J ,t)dm i dai = ^{mi,ai,t)dmidai I dvi -f(r,Vi) / dv Tcl v vcl 'Ei J (v Icl ). (3) 

Jo v i J\ Vi - Vj \ 

The number of collisions between two stars with masses and radii in the ranges (mi, ai) — » (mi + dmi, ai + dai) and (m 3 , aj) — * 
(mj + dm 3 , a 3 + da 3 ) per unit volume per unit time, lZij(r) dm,idaidm 3 da 3 , is given by: 

%3 ( r ,t)dm l da ! dm 3 da 3 ^(^, a ^)dm ! d a! dm 3 d a3 /d 3 v 3 f(r, v 3 )F, (r, v, t), (4) 

and we have TZij = TZji- We define the collision timescale of the stars with mass and radius (mi,ai) as follows: 
. . _ £(mj,ai,t)n(r,t) 

t co n,i (r, t) = , (5) 

J dm 3 da 3 TLji(r,t) 
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where n(r,t) = j d 3 v /(r, v) is the stellar number density at position r and time t. Thus the total number density of the 
stars with mass and radius in the range (mj,Oi) — ► (rrii + drrii,ai + da,i) per unit time undergoing collisions is given by 
n(r, t)£(mi,ai,£)dmidaj/t co ii,i(r, t). For the whole stellar population, the collision timescale t C oii(r) is defined by: 

t ^ r q _ n{r,t) _ 1 ^ 

J drriidaidmjdaj 1Zij(r,t) J dm-idai £(m», Oj, i)/t C oii,i(r,i) 

and the total number of stars per unit time undergoing collisions is given by: 

Note that stellar collision rates depend on the position r in the stellar system and hence vary during the stellar orbit 
(unless the orbit is circular). 



3 FEEDING CENTRAL MASSIVE BHS BY STELLAR COLLISIONS 

In this section, we will study whether the mass released by stellar collisions makes a significant contribution to the growth of 
massive BHs in galactic centers. 

3.1 Stellar mass involved in collisions 

The mass of gas released by two colliding stars depends on their relative velocity and impact parameter, as well as their 
masses and radii and their evolutionary phases (e.g., main-sequence stage or post-main-sequence stage) (e.g., Murphy, Cohn 
& Durisen 1991; Bailey & Davies 1999; Freitag & Benz 2002). We define Amjj(r) as the average gas mass released per 
collision between two stars with masses and radii (mi,ai) and (rrij,aj) at position r. The total gas mass released by stellar 
collisions per unit volume per unit time is given by: 

p co ii(r,t) = i / drriidaidmjdaj TZij(r,t)Amij(r), (8) 



where the collision rate TZij(r, t) can be obtained from equation (4) and we have the factor "1/2" in front of the integration 
because the gas released from each colliding star is counted twice in the integration. The total collisionally released mass per 
unit volume at position r until time t is given by p C oii(r, t) = J Q * p C oii(r, f) dt' , and the total collisionally released mass over 
the age of the stellar system T age is given by: 



M co n(T age ) = / d 3 r p co ii(r,T agc ). (9) 

If all the released gas is accreted onto the central BH, equation (9) gives the mass growth of the central BH caused by stellar 
collisions. If Armj = rm + rrij, equation (9) gives the stellar mass involved or maximum gas mass released in collisions, 

McXCTage). 

To estimate an upper limit on the BH mass growth caused by stellar collisions, we will assume below that all the mass 
involved in collisions is released as gas and accreted onto the central BH. According to equation (6), the variation of the stellar 
number density is given by 

h(r,t) = -n(r,t)/t m n{T,t). (10) 

In equation (10), we have assumed that the stellar density at a position r is affected only by collisions occurring near r, and 
we will see in § 3.3 that our conclusions will not be significantly affected after relaxing this assumption. Now, we assume that 
the stars in the stellar system have a single mass ra» and a single radius a* [i.e., £(m,a,i) = 8(m — m,)S(a — a,), where 6(x) is 
the Dirac function] and the stellar system initially has an isothermal distribution /(r, v) = n(r, 0) exp(— |v| 2 /2cr 2 )/(27rcr 2 ) 3/ ' 2 
(where a is the one-dimensional velocity dispersion) . The initial collision timescale in this stellar system, which can be derived 
from equations (l)-(6), is given by (see also eq. 20 in Duncan & Shapiro 1983): 

*coii(r,0) = — — 2 1 — 2 77-, 2 , r (11) 
16 v / 7rn(r, 0)aia[l + vi BC / {4a 2 )} 

In such a stellar system, if we ignore star formation, stellar evolution, gradual mass segregation or other evolutionary effects, 
then according to equation (6) [i.e., t C oii(r, t) oc l/n(r, t)] and equation (10), we find the stellar number density at time t given 
by: 

""■«>- i+SSbr (I2 > 
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and the stellar collision timescale is given by t co n(r, t) = t + t co ii(r, 0). Thus, in the region where the initial collision timescale 
is much longer than the age of the stellar system [i.e., t co n(r,0) 3> Tagc], we have n(r, T agc ) ~ n(r, 0) (i.e., the stellar 
number density is affected little by collisions) and the collision timescale i co ii(r, T ago ) — i C oii(r,0); while, in the region with 
icoii(r,0) <C Tkgo, we have n(r, T agc ) ^ n(r, 0)i C oii(r, 0)/T agc <C n(r, 0) (i.e., stars are disrupted quickly by stellar collisions) 
and the collision timescale is t co ii(r, T ago ) — ^age- The total gas mass released by collisions (or the upper limit of the total 
stellar mass involved in collisions) per unit volume over the age of the stellar system is given by: 

P™(r^ = T agc ) = m„[n(r,0) - n(r,T age )] 

m»7i(r,0)T a goAcoii(r, 0) if t co ii(r, 0) > T agc 

m*n(r,0) if i co ii(r, 0) < T ag e ' 1 ' 

and the total mass involved in stellar collisions over time T age is given by: 

M c X(Tage) = j dr 3 Pc ^r(r,T ago ) (14) 
If the initial stellar density is singular, that is, 

n ^=2^ < 15 > 

(cf., eq. 4-123 in Binney & Tremaine 1987), then using equations (11)— (15), we have the total mass involved in stellar collisions 
over time T age as follows: 

M™?r(T ag e) = VAt'/G, (16) 

where 



A = GfZ \ 1+ 4^)- < 17 > 

In a stellar system with a distribution of stellar masses and radii, if we ignore star formation, stellar evolution, gradual 
mass segregation or other evolutionary effects and we assume that collision timescales of stars with any given mass and radius 
(mj,Oj) do not change with time, i.e., t C oii,i(f, t) = t co n,i(r, 0), we may obtain the number density of stars with mass and 
radius (mj,Oj) at time t is H(m.j, 0)n(r, 0)e _t/ ' tcoll ' i ' r ' - ) (cf., eq. 10), and the stellar mass involved in collisions over the age of 
the stellar system is given as follows: 



ACjT(Tage) = J 47rr 2 dr J dm, m l H(m,, 0)n(r, 0) [l - e - T **-/*«*.<W] . (18) 

In reality, the collision timescale £ C oii,i( r ) t) is n °t a constant and should increase with time, M™ 1 a 1 x (T agc ) obtained from equation 
(18) is an upper limit to the stellar mass involved in collisions. 

Ignoring star formation and stellar evolution in the study is usually a reasonable assumption for an old stellar system 
after the initial rapid stellar evolution of ~ 1 Gyr is complete (i.e., massive stars have evolved into stellar remnants and the 
main-sequence stars have low masses with long lifetime and little change in the mass and radius in the main-sequence phase; 
cf., § 3.4). Stellar remnants of massive stars (BHs and neutron stars) may migrate inwards as a result of energy equipartition 
with low-mass field stars (e.g., Morris 1993), i.e., mass segregation. As will be discussed at the end of § 3.4, mass segregation 
affects little on our result of M^f* in most galactic centers, especially for bright galaxies. 



3.2 Galaxy samples 

In the past several years, images from the Hubble Space Telescope (HST) have revealed many details about the central regions 
of nearby galaxies, with a resolution of 0.1", corresponding to distances of < lOpc or ~ 10 6 (M./10 8 Mq) Schwarzschild radii 
for typical target galaxies (Byun et al. 1996). The inner surface brightness profiles of the galaxies as a function of projected 
radius 7? are well fitted with a five-parameter fitting function — the Nuker law (Faber et al. 1997): 



m =2^i b (-y [i + ( 



(19) 



The asymptotic logarithmic slope inside r b is —7, the asymptotic outer slope is — /?, and the parameter a characterizes the 
sharpness of the break. The break radius r b is the point of maximum curvature in log-log coordinates. The "break surface 
brightness" I b is the surface brightness at r b - According to the inner surface brightness profiles I(R) oc RT 1 (r — » 0), elliptical 
and spiral bulges (hot galaxies) can be classified into two types: core galaxies (7 < 0.3) and power-law galaxies (7 > 0.5). 
A sample of nearby hot galaxies (elliptical galaxies or spiral bulges) observed by HST is collected in Faber et al. (1997), 
who estimate the stellar mass-to-light ratio T (constant for each galaxy, determined by normalizing to the central velocity 
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dispersion based on spherical and isotropic models fitted to the Nuker-law profile) and the half-light radius (cf., Table f in 
Yu 2002). Assuming that the galaxies are spherical and isotropic (the mean ellipticity is 0.26), we may obtain their stellar 
density and velocity distributions by using the Eddington formula (eq. 4-f40a in Binney & Tremaine 1987) and estimate 
their central BH masses M, by using the following relation between BH mass and galactic velocity dispersion (Tremaine et 
al. 2002) 

log(M./M ) = (8.13 ± 0.06) + (4.02 ± 0.32) log(a e /200 kms -1 ), (20) 

where u c is the luminosity-weighted line-of-sight velocity dispersion inside the half-light radius. Note that the surface brightness 
profiles obtained in Faber et al. (1997) are based on HST WFPC1 data. We have updated the surface brightness profile of 
M32 using WFPC2 data in this paper (for detailed information, see § 4.2 or Lauer et al. 1998). In total, 31 galaxies from the 
Faber et al. sample will be used in this paper to study the effects of stellar collisions in galactic centers. 



3.3 Results for simple stellar systems with identical stars 

In this subsection, we use equation (6) to find the collision timescales as a function of galactic radius in realistic galactic 
centers (cf., § 3.2). The results are shown in Figure 1. Then, we will use equations (8)-(17) to obtain the stellar mass involved 
in collisions or the upper limit of mass released by collisions (M™^) over the Hubble time tHubbio — 10 10 yr, which is shown 
in Figure 3. In this subsection, all the stars are assumed to have identical mass and radius and stellar evolution is ignored. In 
the next subsection, we will see that the generalization to a distribution of masses and radii will not significantly affect our 
conclusions. 

Figure 1 shows the collision timescales i C oii,©(r) as a function of radius r (the subscript "0" indicates that the results are 
obtained by assuming that all stars have the solar mass and radius). The stellar distribution n(r) and the collision timescale 
tcoii,©(f) within the region not resolved by the HST (2r < 0.1") are obtained by extrapolating the observed surface brightness 
profile (eq. 19) inward [i.e., assuming I(R) oc RT" 1 at 2R < 0.1"]. As seen from Figure 1, the collision timescales t C oii,©(f) 
increase with increasing radius r. At radii r > 0.003-1 pc (Fig. la), the collision timescales are longer than the Hubble time 
and the stellar distributions are affected little by stellar collisions [i.e., n(r, 0) ~ n(r, teubbie) and i co ii(r, 0) ~ t C oii(r, teubbie)], 
which is also true in the region resolved by HST (2r > 0.1", cf., Fig. lb). Figure 1(c) shows that in all galaxies, the region 
with £ C oii,oM < iHubbie is located within the radius of the sphere of influence of the BH m, which is defined in terms of the 
intrinsic one-dimensional velocity dispersions of the galaxy a(r) through 

a 2 {r H ) = GM./r H . (21) 

In the region where most stars should have been disrupted quickly [i.e., t C oii,©(f) "C iHubbie], the stellar distribution obtained 
by extrapolating the Nuker law inward can only be interpreted as the initial stellar distribution and the corresponding collision 
timescales t C oii,©(f) shown in Figure 1 are the initial collision timescales t co ii,©(r, 0). Current stellar densities in those regions 
may be estimated from equation (12), and current collision timescales are about t co n(r, tHubbie) — ^Hubble- 

Figure 2 shows the collision timescale £ C oii,© as a function of the ratio 2a(r)/v esc ,o, where iw,© is the escape velocity 
from stars with solar mass and radius (cf., eq. 1). As seen from Figure 2, the one-dimensional velocity dispersion a(r) is not 
a monotonic function of i C oii,©(r) or r (note that t co ii,© is a monotonic increasing function of r in Fig. 1), which is due to 
the effect of central BHs. In the absence of central BHs, the velocity dispersion a(r) should decrease with decreasing radius 
r at small r (cf., eq. 10 for case 1 < ij < 2 in Tremaine et al. 1994); but with a central BH, the velocity dispersion increases 
with decreasing r at r C th [i.e. cr(r) oc r -1 / 2 ]. Thus, with decreasing r or t co \\ t ©(r), the velocity dispersion a(r) may show 
a minimum near the region where the effect of central BHs becomes dominant. Figure 2 shows that for all the core galaxies, 
the collision time is less than the Hubble time only if the velocity dispersion exceeds the escape speed. 

Figure 3(a) shows the central BH mass M. versus the upper limit of the mass released by stellar collisions over a Hubble 
time JVf "ff^Q^Hubbie) (eq. 9) in the galaxies in our sample (open circles and solid circles). Once again, all stars are assumed to 
have the solar mass and radius. As seen from Figure 3(a), the upper limit of the mass released by stellar collisions over a Hubble 
time M^f^Q (inubbic) is smaller than the BH mass M, in all the galaxies of the sample. For power-law galaxies (M. ~ 10 5 - 
1O 9 M0, open circles), M™f* is usually in the range ~ 1O 4 -1O 7 M0; and for core galaxies (M. ~ 1O 8 -1O 9 M0, solid circles), 
M™^ Q ~ 10 4 -10 6 M Q . Thus, other material sources, rather than mass released by stellar collisions, must dominate the growth 
of central BHs, especially in core galaxies. In Figure 3(a), we also show the results obtained from simplified galaxy models with 
Maxwellian stellar distributions with dispersion independent of radius (dashed lines). The short dashed lines are obtained by 
assuming that the galaxies have flat inner stellar densities (see eqs. 4-124b and 4-128a in Binney & Tremaine 1987) 

n(r) = n c /[l + (r/r c ) 2 ] 3/2 , (22) 

where n c = 9<r 2 /(47rGm»r 2 ). The core radius r c is set to be in the range 10 2 -10 3 pc. Note that it is only in the singular case 
that the density in equation (15) and a Maxwellian velocity distribution are an exact solution of the collisionless Boltzmann 
equation. The long dashed line is obtained by assuming that the galaxies have singular stellar densities (see eqs. 15-17). 
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Figure 1. Stellar collision timescales as a function of galactic radius (eq. 6). The galaxy sample comes from Table 1 in Yu (2002). Stars 
are assumed to all have the solar mass and radius. In the region not resolved by the HST (2r < 0.1"), the collision timescales are obtained 
by extrapolating the surface brightness profiles (eq. 19) inward. The solid curves represent core galaxies (7 < 0.3) and the dotted curves 
represent power-law galaxies (7 > 0.5). The thick dashed line represents the galaxy M32, whose central colors will be studied in detail 
in § 4.2. The three panels use three different measures of radius: parsecs, arcsec and rjj, where rjj is the radius of the sphere of influence 
of the BH, defined implicitly in terms of the intrinsic one-dimensional velocity dispersion of the galaxy <r(r) through equation (21). For 
how the results shown in this Figure are affected after generalizing to a distribution of masses and radii, see § 3.4. 



© 0000 RAS, MNRAS 000, 000-000 



Stellar collisions in galactic centers: black hole growth and color gradients 7 




Figure 2. Stellar collision timescales t C oll,0 as a function of the ratio cr(r)/v eBCi Q. The quantity cr(r) is the one-dimensional velocity 
dispersion and « e sc,© — 620 km s -1 is the escape velocity from stars with solar mass and radius (cf., eq. 1). The galaxies are the same as 
those in Figure 1. Stars are assumed to all have the solar mass and radius. The line types have the same meaning as those in Figure 1. 
This figure shows that the velocity dispersion <r(r) is not a monotonic function of radius r, which is due to the effect of central BHs (note 
that t coU ^Q is a monotonic function of r in Fig. 1). 



Here, the isothermal distributions with the flat and singular stellar densities represent simplified cases for core and power-law 
galaxies, respectively. As seen from Figure 3(a), most of the solid circles (obtained from actual core galaxies) are located in 
the region covered by the short dashed lines (those solid circles not covered have break radii rb out of the assumed range 
of r c ~ 10 2 -10 3 pc). The mass involved in collisions obtained from the singular isothermal distribution (long dashed line), is 
basically an upper bound to the masses obtained from actual power-law galaxies (open circles), because the singular stellar 
density profile is generally steeper than the inner stellar density profiles of most power-law galaxies. 
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Figure 3. Panel (a) shows the total stellar mass involved in collisions M™^ X 0(*Hubble) ( e q- 9) over a Hubble time, versus the central 
BH mass M, . The galaxies are the same as those in Figure 1. The BH mass M. is obtained by the BH mass versus galactic velocity 
dispersion relation (eq. 20) in Tremaine et al. (2002). Stars are assumed to have the solar mass and radius. The stellar mass involved in 
collisions is an upper limit to the mass that collisions can contribute to the BH growth (assuming both stars are completely disrupted and 
the entire gas mass is accreted by the BH). The solid circles represent core galaxies; and the open circles represent power-law galaxies. 
The dotted lines arc the reference lines representing M™^ X (*Hubblc) = Mm and M™J*[ X (tHubblc) = 0.1M.. The dashed curves give 
the relation between M. and M™J t 1 x Q (tH u bble) m simplified galaxy models with isothermal stellar distributions. The long dashed curve 
gives the results with singular stellar densities (see eqs. 15-17). The short dashed curves give the results with flat inner stellar densities 
(eq. 22); reading upwards, the core radius r c in equation (22) varies from 10 2 pc to 10 3 pc with interval Alog(r c /pc) = 0.2. In all cases 
M™j l 1 x Q (tHubblc) < Mm. Panel (b) shows 47rr 3 p™^ X g (r, tHubblc) as a function of radius r (cf., cqs. 8 and 9). The peaks of the curves give 
the location where most of the stellar mass involved in collisions originates. The line types have the same meaning as those in Figure 1. 
For how the results shown in this Figure change after generalizing to a distribution of masses and radii, see § 3.4. 
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Figure 3(b) shows 47rr 3 p™n^ (r, t Hu bbio) (cf., eqs. 8 and 9) as a function of galactic radius r. The peaks or highest parts 
of these curves represent the regions where most of the stellar mass involved in collisions comes from. The radii where the 
peaks of the curves are located can be denoted by r pca k. We have (v 2 e i) = 6er 2 (r) for the typical relative velocity of two 
colliding stars at radius r. As seen from Figure 3(b), for core galaxies (solid lines), the peaks of 47rr 3 p™^ (r) are generally 
located at large radii (r poa k ~ 10 2 -10 3 pc). For power-law galaxies (dotted lines), some of the curves show apparent peaks at 
small radii (r pca k ~ 0.1-1 pc); while some have flat-topped profiles covering a large range of radii (e.g., 10~ 2 -10 2 pc). Note 
that the mass-to-light ratios of the galaxies are obtained by normalizing to the central velocity dispersion based on spherical 
and isotropic models fitted to the Nuker-law surface brightness profiles. Thus, given the surface brightness profile and the 
mass-to-light ratio of the stellar system, we have at r < rb (tv the break radius, cf., eq. 19) n(r) oc r _1 /(r)/m* oc r _7_1 /m* 
for power-law galaxies. Now, consider two regions with r <C rb in power-law galaxies: 

• In a region with t co n(r) < tHubbie, we have 

r 3 p™n x (r) oc r s n(r)m* oc r~ 7+2 , (23) 

and since 7 < 2, r 3 p™^ x (r) increases with increasing galactic radii r (this applies to the region r < 0.01-0.1 pc in Figure 3b). 

• In a region with t C oii(r) 3> tHubbie, we have (from eqs. 1, 6 and 19) 

. , , / r y+1 a(r)~ 1 m*/al if cr(r) > t> csc . . 

t C oii(r) oc [n r )a r oc i +1 ) ' )' 24 

[ r 1 ^ a(r)/a* if a(r) < v csc , 

and 

'*™ « ^/u W «{;:::;j;{*' ;-;:>>> ::: m 

The (one-dimensional) velocity dispersion cr(r) varies only slowly with the galactic radius r at r > ru, and varies as r -1 ^ 2 
at r < r H , where r H is the radius of the sphere of influence of the BH (cf., eq. 21). In the region with er(r) 3> v eBC (r), we 
usually have r <C rH- Thus, according to equation (25), with increasing r, r 3 p™"ff(r) monotonically decreases with radius for 
power-law galaxies with 7 > 0.75. For power-law galaxies with 0.75 > 7 > 0.5, r 3 p™f x (r) is usually not a monotonic increasing 
or decreasing function of r due to the different variation of a(r) with r at r < ru and r > ru- 

For power-law galaxies with 7 > 0.75, since r p™^ x (r) is an increasing function of r in the region with t C oii(r) <C tHubbie and is a 
decreasing function of r in the region r < r\, with t co n{r) ^> tHubbie, r 3 Pcoif ( r ) shows a peak in the region with t co u(r ) ^ teubbie, 
and the mass involved in collisions comes mainly from small radii r pea k ~0.1-lpc (cf., Fig. la) (see Figure 3b). However, 
for core galaxies, the peak of r 3 p™^ x (r) is generally not located around the radii with t co n(r) ~ tHubbie, but at large radii 
r P eak ~ 10 2 -10 3 pc where the stellar densities show a break (see the break radii r b in eq. 19). In our analysis below, we 
will focus on power-law galaxies with 7 > 0.75 and core galaxies, which usually can be analytically studied; the results for 
power-law galaxies with 0.75 > 7 > 0.5 are usually intermediate cases between the results of those two types of galaxies. 

Note that in equation (10), we assume that the stellar density n(r) is affected only by local stellar collisions. This 
assumption will not underestimate the mass involved in collisions per unit volume at position r if the collision timescale at 
position r is much longer than the age of the stellar system, and this assumption will not affect the estimate of the mass involved 
in collisions per unit volume at position r if stars are on circular orbits. For core galaxies, the mass involved in collisions 
M™f x comes mainly from near the break radii rb, where t C oii(r b ,0) ^> tHubbie; hence, M™f x ~ m*n(rb)tHubbie/tcoii(rb, 0) will 
not be affected even if we relax the assumption in equation (10). For power-law galaxies with 7 > 0.75, the mass involved 
in collisions M™^ comes mainly from the radii r pcak where t C oii(r pea k) — tHubbie, and M™fi x is approximately equal to the 
initial stellar mass within the radius r pca k (cf., eqs. 13 and 14). Note that the stars on eccentric orbits which are not involved 
in collisions at their apocenters (e.g., > r pea k), may be involved in collisions at their pericenters (e.g., < r pea k), where the 
collision timescale is shorter due to high stellar density; and the mass involved in collisions in Figure 3(a) is therefore possibly 
an underestimate. However, the total stellar mass involved in collisions at radius r cannot exceed the total mass of stars 
which may come within radius r. To find an upper limit to the correction needed for non-locality, we assume that all the 
stars are involved in collisions at their pericenters (if the collision timescale at their pericenter is t co n(r) < tHubbie) and the 
total mass involved in collisions per unit volume at position r is about min[m»n(r, 0)tH ub bie/tcoii(f, 0), (l/47rr 2 )dM< r (r)/dr], 
where M< r (r) represents the total mass of stars which may come within radius r in spherical stellar systems (cf., Fig. la in 
Yu 2002). We find that with these assumptions, M™f x for power-law galaxies with 7 > 0.75 in Figure 3(a) will increase at 
most by a factor of 2-5 (average 3). 

The collision timescale and the stellar mass involved in collisions shown in Figures 1 and 3 are obtained by assuming that 
all of the stars have the solar mass and radius. Now consider the more general case in which all the stars have identical mass 
and radius (m„, a*), which are not necessarily the solar values. We find the following changes in the results shown in Figures 
1 and 3: 
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0.1 M 


, 0.16R© 


10 M 


3 , 6.3 R Q 




core galaxy 


power-law galaxy 
(0.75 <7 < 1.0) 


core galaxy 


power-law galaxy 
(0.75 <7 < 1.0) 


^peak/^pcak ,0 


1 

0.2 


0.2-0.6 
0.2-0.6 


1 
6 


1.7-4.4 
1.7-6.3 


tcollAcoll,© 




4-6 





2-0.3 



Table 1. Changes in the results shown in Figures 1 and 3 (the stellar collision timescale t co ii, the total stellar mass involved in collisions 
M™j5 x , and the radius where most of the stellar mass involved in collisions originates r pea k) by assuming that all the stars have mass 
and radius lower than the solar values (0.1 Mq, 0.16 Rq) or higher than the solar values (10 Mq, 6.3 Rq). See § 3.3. 

• The collision timescales vary as (m*/ Mq)( Rq/o,) 2 [if v eBC <JC a(r)] or Ro/a« [if v CBC ^> a{r)\ times i C oii,©(r) shown in 
Figure 1 (cf., eq. 24). 

• In the region with t co n(r) <JC ^Hubble, 47rr 3 p™^ x (r) is the same as that shown in Figure 3(b) (cf., eq. 23). In the region with 
tcoii(r) > teubbic, 47rr 3 p™if(r) varies as (o»/ R ) 2 ( M Q /m») [if v esc < a(r)] or a*/Rg [if v csc 3> <r(f)] times 4irr p^,n^©( r ) 
shown in Figure 3(b) (cf., eq. 25). 

• For core galaxies, the mass involved in collisions M™f 1 x (tHubbi e ), which comes mainly from the region near the break 
radius rt,, should be gs*/Rq times that shown in Figure 3(a). 

• For power- law galaxies with 7 > 0.75, the mass involved in collisions M^j^iHubbic) comes mainly from the region 
r ~ r poak < r H with icon(r pe ak) — taubbic; thus, according to equation (24), r pcak should be [(o»/ R©) 2 ( M Q /m*)] 1/(7+1B) 
[if a(r) > Uosc] or (a„/R Q ) 1/(7+0 - 5) [if a(r) < v csc ] times that shown in Figure 3(b), and M™g* [oc r 3 cak p™ti x (r pcak ) oc 
^ cak n(r pcak )/Wbic oc r p 7+ 2 ] is about [(a*/ R Q ) 2 ( Me/m,)]^ 2 ^ 1 ' 5 ) [if a(r) » v csc ] or (a./ R )(- 7 + 2 >/( 7 +O-5) [if a (r) « 
v csc ] times M™" shown in Figure 3(a). 

We now give two examples of the above changes. We first assume that the galaxies are composed of low-mass stars with 
m* = 0.1 M© and a, = O.16R [~ (m./Mg)' 7 , 77 = 0.8, Kippenhahn & Weigert 1990]. The collision timescales t co n will 
increase by a factor of ~ 4 [v CBC <JC o(r)] or ~ 6 [v csc >• a(r)] compared to those in Figure 1. For core galaxies, the radii 
where the peaks of 47tr 3 p™^ x (r) are located, r poak , are generally still around the break radii r b , but the stellar mass involved 
in collisions M™f x will decrease by a factor of ~ 0.2 compared to Figure 3. For power-law galaxies with 1.0 > 7 > 0.75, r pcak 
will decrease by a factor of ~ 0.2-0.6, and so will the stellar masses involved in collisions M™f x . Similarly, we may also obtain 
the changes by assuming that the galaxies are composed of high-mass stars with m* = 10 Mq and a* = 6.3 R . The results 
are summarized in Table 1. 



3.4 Generalization to a distribution of stellar masses and radii 

A realistic stellar system is composed of stars with a distribution of masses and radii, rather than identical stars. To investigate 
this case, we assume that all stars are formed instantaneously at the formation of the stellar system, and use the following 
two forms of stellar initial mass function (IMF): one is the Salpeter IMF (Salpeter 1955), 



n u , _ - / ^sp(m/ M )- 2 - 35 0.08 < ml M Q < 120, 
otherwise; 

and the other is the multi-power-law IMF (cf., Kroupa 2002), 

(m/O.O8M )- ' 3 0.01 ^ m/ Mq < 0.08, 



(26) 



H mu iti(m, t — 0) — A m uiti < 



(m/O.OSM©)- 1 ' 3 0.08 ^ m/Mg < 0.5, (2? n 

(O.5M Q /O.O8M )- L3 (m/O.5M )- 13 0.5 < m/ M Q < 120.0, 
otherwise. 

The constants As p and A mu i ti in equations (26) and (27) are determined by J °° H(m, t = 0) dm = 1. The multi-power-law 
IMF gives a much flatter slope at the low-mass end than the Salpeter IMF. For stars with initial mass m, the lifetime on the 
main sequence is about (eq. 9-3 in Binney & Tremaine 1987) 

t MS ^ 10Gyr(m/M Q )- 2 ' 5 . (28) 

In a stellar system with age T ago ~ 10 Gyr, the mass of stars at the turn-off point m t f is about IMq. The present mass 
function of the remaining low-mass main-sequence stars follows their initial mass function H(m, 0) with m < m t f . These stars 
are assumed to have radii a ~ (m/MQ^Rg with r\ = 0.8 (p. 208 in Kippenhahn & Weigert 1990). Massive stars (> IMq) 
usually have lost a significant fraction of their mass and become low-mass remnants (BHs, neutron stars or white dwarfs). We 
assume that progenitors with mass m > 30 Mq will become BHs with mass 8 Mq , progenitors with mass 30 Mq > m > 8 M 
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will become neutron stars with mass 1.5 M and radius 10 km, and progenitors with mass 8M > m > 1M Q will become 
white dwarfs with mass 0.6Mq and radius 10 4 km (e.g., Murphy, Cohn & Durisen 1991). Besides stellar remnants, there also 
exist other, relatively short lived, post-main-sequence stars (such as red giants, horizontal branch stars or asymptotic giant 
branch stars etc.) which have evolved from stars with initial mass close to mtf. 

According to equation (5), we have the collision timescale of main-sequence stars with mass and radius (m», a t ) as follows: 

tcoii,,(r,r agc )oc 1 . (29) 

n(r,T age ) J dmjdaj (,{mj,aj,Ta %e )Y, ii 

Given the dynamically determined mass-to-light ratio (by normalizing to the central velocity dispersion based on spherical 
and isotropic models fitted to the surface-brightness profile, see § 3.2) and the observed surface brightness profile, we have the 
following relation for the stellar number density obtained by using the Eddington formula (eq. 4- 140a in Binney & Tremaine 
1987) 

n(r,T ago ) oc l — — — , (30) 

J () dmj m 7 -a(mj,Tage) 

which depends on the stellar mass function H(mj,T ago ). Thus, using equations (1), (29) and (30), we have 

^coll,i(^*, ^k g c) 
£coll,© (?*] Tage) 

dmj (mj/M )H(mj,T agc ) 



drrij H(m J ,O)[(a l + a i ) 2 /(2R ) 2 ] if a(r) > u esCiQ 

dm, H(mj , 0) [(% + a 3 )/2 R ] [(m, + mj)/2 M ] if a(r) « u esCiQ ' 1 ' 



where t co \\,o (see § 3.3 or Fig. 1) is the collision timescale obtained by assuming that all the stars have the solar mass and 
radius, and t> C sc,© = y/2G M© / R© is the escape velocity from stars with solar mass and radius. The two cases in equation 
(31) are divided by v eS c,e ~ ff(r) (cf., eq. 1), because we have Vij = y2G(m, + mjj/jai + Oj) ~ v eS c,o for any collisions 
of two main-sequence stars with masses m< and rrij in the range 0.01-1 M© [assuming a oc mP with rr(=0.8) being close to 
1] and we assume that the mean relative velocity of two colliding stars is \J (v TC \) (v~^)~ x ~ ( v t C i) 1 ^ 2 = \f§o(r) ~ <r(r) (see 
eq. 1). In equation (29), we have ignored collisions with stellar remnants (BHs, neutron stars or white dwarfs), because the 
number of stellar remnants (evolved from massive stars) is small compared to that of main-sequence stars (note that mass 
segregation is ignored) and also the radius of stellar remnants is not large (< 1R©). We have also ignored collisions with 
other post-main-sequence stars (e.g., red giants) since the number of these stars is also small. Note that the large radius of 
giants (e.g., ~ 10 2 R ) will increase the collision rates of main-sequence stars, but collisions with red-giant envelopes generally 
result in little mass loss from either star because the amount of mass contained in the greatly extended envelopes of giants is 
typically a small fraction of a solar mass and the overlapping area of collisions is quite small relative to the total collision cross 
section. In addition, as argued by Murphy, Cohn & Durisen (1991), giants undergo significant mass loss in a short amount 
of time anyway because of stellar winds, and it makes not much difference whether the mass comes from winds or collisions 
in the sense of providing raw material for the BH growth. 

Using equation (29) and assuming that the stellar radius m = (rm/ M©)~ 0,8 R© and the mass of stars at the turn-off 
point mtf = IMq, we obtain the collision timescales t co u,i as a function of stellar mass rm. The ratio of collision timescales 
£coii,i/£coii,0 (eq. 31) as a function of stellar mass rrn is shown in Figure 4 (see t C oii,© in Fig. 1). We show the results obtained 
both by assuming the Salpeter IMF (eq. 26) and by assuming the multi-power-law IMF (eq. 27) . As seen from Figure 4, the 
ratio icoii.i Acoii.© increases from ~ 0.6 to ~ 7 (or to ~ 9) with decreasing stellar mass rrii, from m t f to 0.08 M for the Salpeter 
IMF (or to 0.01 M for the multi-power-law IMF), and the difference of the ratios between the two IMFs for a(r) <JC v eS c,@ 
and for a(r) 3> f CS c,© is less than 30 percent. For the whole population of main-sequence stars, we use equation (6) and the 
collision timescale t C oii,i shown in Figure 4 to obtain the total collision timescale t co ii, which is about 3-4 times the timescales 
tcoii,© shown in Figure 1. 

According to equation (8), we have the total mass density of main-sequence stars involved in collisions with main-sequence 
stars with the natural logarithm of their masses in the range In m.j — > In m, + d In m; per unit time given by: 

r m tf 

PraiM( r . T a gc )dlnm 4 oc dm* / dm^- (m, + m i )7l, J (r,r ilgc ) 

Jo 
r m tf 

dm, / drrij (m-i + mj)H(mi,0)H(m,j,0)n 2 (r,T age )Eij; (32) 



oc 



o 



and using the similar derivation to obtain equation (31), we have the ratio 



Pcoll,i( r i Tag< 



miE(rm, 0) 



dmj (mj/M )S(mj,T E 



a g cj 



-2 
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/""dm, S(m j ,O)[(a, + o ; )/(2R )] 2 [(m, + ra j )/(2M a )] a(r) » v, 



esc,© 



(33) 



/""dmj H(m,,0)[( ai +a j )/(2R )][(m i + m J )/(2 M Q )] 2 a(r) < tw,© 

where p™n x © (cf., eq. 8) is the total mass density of main-sequence stars involved in collisions per unit time obtained 
by assuming that all the stars in the stellar system have solar mass and radius. We put the factor "1/2" in front of 
pS>im('"> ^a g c)/p™n^0( r ) lago) because / Q mtf P™?u( r > ^age) dlnm; gives twice of the mass involved in collisions (see also the 
factor "1/2" in eq. 8). We also have the total mass density of main-sequence stars with the natural logarithm of their masses 
in the range In rrii — > In rrii + d In rrii involved in collisions per unit time: 

PcoTi?r(r,ra gc )dlnmi oc dm, / dm, rrnTZtjir,^^) 

Jo 

oc dnii / dm, miH(m», 0)E(m 3 ;, 0)n 2 (r, T ago )Eij, (34) 



and 

■ / max / rp 
Pcoll.i V' ! J a ( 



(m 4 /M )m l H(m i ,O) 



dm, (mj/M Q )H(m J ',r £ 



(35) 



Proll*© ( r i^agc) 

J mtt dm,- H(m,-,0)[(oi + a,-)/(2R )] 2 cr(r) > tw,o 

Jo dm,- S(mj,0)[(oi +a j )/(2R )][(m, +m,)/(2M )] <j(r) < iw,© 

As in equations (29) and (31), collisions with stellar remnants and other post-main-sequence stars are ignored in equations 
(32)-(35). We use equations (32)-(35) to obtain p^Si ti (r,T agc ) and p' c ™{^?(r, T ag c) as a function of stellar mass m». The ratios 
of p^u/(2p m >n X 0) and p^cST /Pcoif© as a function of stellar mass m, for both the cases o(r) <C ?w,q (dotted lines) and 
<r(r) >• Wcsc.q (solid lines) are shown in Figure 5(a) and (b), respectively. We also show the results obtained both by assuming 
the Salpeter IMF (eq. 26) and by assuming the multi-power-law IMF (eq. 27). In Figure 5(a), the curves of p™nV(2p™©) 
obtained from the Salpeter IMF (with the low-mass end at 0.08 M©) show a minimum at an intermediate mass between 
0.08-1 M© and increase from the minima to both the the high-mass end (~ 1M©) and the low-mass end (~ 0.08 M©), which 
reflects the fact that the stellar mass involved in collisions comes mainly from collisions between high-mass stars and low- 
mass stars; and the curves of p™n X i/(2pSS,n x ©) obtained from the multi-power-law IMF (with the low-mass end at 0.01 M©) 
increase with increasing stellar mass rrii, which reflects the fact that the stellar mass involved in collisions comes mainly from 
collisions with high-mass stars since the multi-power-law IMF gives a natter slope at the low-mass end than the Salpeter 
IMF. Figure 5(b) shows that PcoU^f I 'Pmii,® obtained from both the Salpeter IMF and the multi-power-law IMF increase with 
increasing stellar mass rrii, which reflects that the stellar mass involved in collisions comes mainly from high-mass stars, 
especially for the multi-power-law IMF since the slope of corresponding curves in Figure 5(b) is steeper. For core galaxies, 
the total mass of main-sequence stars involved in collisions comes mainly from galactic radii r ~ [where t cc ,ii,i(r) 3> ^Hubble 
and cr(r) <JC t> osc ,©], which is found to be ~ f™ t! P™oi\,i d lnmi/(2p™^ x ©) ~ 0.4-0.5 times the result obtained by assuming that 
the stellar systems have identical stars with solar mass and radius. For power-law galaxies, using equation (18), we find that 
the total stellar mass involved in collisions M™f x is smaller than M™j5 x © (see § 3.3 or Fig. 3) by a factor of 0.4-0.7. 

As shown above, after the generalization to a distribution of stellar masses and radii, the total stellar masses involved in 
collisions decrease compared to those shown in Figure 3(a). In addition, note that not all the collisions result in disruption of 
stars and not all the gas released by collisions can be accreted onto BHs. Hence, the conclusion obtained from Figure 3 that 
the mass released by stellar collisions is not enough to account for the central BH growth (especially in bright core galaxies) 
will not be changed, by more realistic models. 

Mass segregation is ignored in the analysis above. This assumption affects little on our result of stellar mass involved in 
collisions at least for all the core galaxies and most of the power-law galaxies for the following reasons. Our calculation shows 
that the relaxation timescales of stellar remnants of massive stars (BHs and neutron stars) are higher than the Hubble time 
for all the core galaxies and are higher than the Hubble time at r > rp Ca k,© for most of the power-law galaxies [only for about 
four power-law galaxies, the relaxation timescales of the stellar BHs (8M©) are smaller than or approximately equal to the 
Hubble time at the region r ~ r peaki ©.]. 



3.5 Material sources for growth of central BHs 

As mentioned in § 1, gaseous material for the growth of central BHs in an isolated stellar system may come from stars through 
three mechanisms (e.g., Murphy, Cohn & Durisen 1991): tidal disruption (or being swallowed whole by central BHs) of stars 
(mostly on elongated radial orbits), stellar collisions, and stellar evolution. In this subsection, we will compare these three 
contributors to BH growth. 

In a stellar system with a central BH, stars that come within a distance r tid ~ (M./m») 1 ' /3 a» of the central BH will be 
tidally disrupted by the BH if r tid > rs c h (where rs c h = 2GM,/c 2 is the Schwarzschild radius) or swallowed whole by the BH 
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if f*tid < fsch- Here, we do not distinguish the above two cases and simply call the total mass of stars which may come within a 
distance of r — max(r t id, 7*s c h) as the tidally disrupted mass M t ^. We use the same method as Magorrian & Tremaine (1999) 
to obtain Mud in realistic galaxies, and the galaxies are assumed to be spherical since flattening or triaxiality of galaxies 
does not significantly change the results (Magorrian & Tremaine 1999). Figure 6 shows the tidally disrupted stellar mass 
(Mud,©) versus the stellar mass involved in collisions (M^^q) over a Hubble time; in this Figure, the stellar systems are 
assumed to be composed of identical stars with solar mass and radius. As seen from Figures 3 and 6, M t i d>0 and M™^ Q are 
generally comparable. After the generalization to a distribution of stellar masses and radii (i.e., the Salpeter IMF in eq. 26 or 
the multi-power-law IMF in eq. 27), M tid will increase by a factor of ~2-3 (cf., Magorrian & Tremaine 1999), while M^f* 
will be reduced by a factor of ~ 0.4-0.7. Considering that not all the collisions result in stellar disruption and not all the gas 
released by collisions can be accreted onto BHs, tidal disruption is likely to contribute more mass to the central BH growth 
than stellar collisions. Magorrian & Tremaine (1999) show that tidal disruption contributes significantly to the present BH 
mass only in faint galaxies < 10 9 L Q . For bright galaxies, our study in this paper shows that neither tidal disruption nor 
stellar collisions can contribute significantly to their present central BH masses. 

Using the simple model of stellar evolution described at the beginning of § 3.4, we may obtain the mass loss caused by 
stellar evolution if all stars are formed instantaneously at the formation of the stellar system. The mass loss may comes either 
from stellar winds (mainly during giant stages) or from the mass lost by Supernovae II. The total mass loss amounts to about 
21 percent of the initial total stellar mass in the initial 1 Gyr and about 7 percent during the interval from 1 Gyr to 10 Gyr for 
the Salpeter IMF. The corresponding numbers for the multi-power-law IMF are 32% and 10%. (For both of the IMFs, more 
than forty percent of the mass loss in the initial 1 Gyr comes from stars with initial mass in the range 2.5-8 M Q or mass lost 
by stellar winds.) Therefore, mass loss during stellar evolution at the early stage of galaxy evolution is a possible source to 
feed central BHs if part of the lost mass can lose angular momentum and sink to galactic nuclei. 

A final possibility for the material source is that galaxy encounters may concentrate gas in the centers of interacting 
galaxies to feed BHs (see simulations in Barnes & Hernquist 1996). The gas infall may fuel pre-existing central BHs; or 
galaxy encounters may trigger a central star burst, which could lead to nuclear activity. For both of the above possibilities, 
massive BHs form or grow rapidly, early in the universe's history (e.g., at redshift z ~2-3 when galaxy encounters mostly occur 
in the hierarchical galaxy formation model). This scenario is consistent with one of the arguments made by Yu & Tremaine 
(2002) that growth of high-mass BHs (> 10 Mq) comes mainly from accretion during optically bright QSO phases (which is 
obtained by studying the relation between the local BH mass function and the QSO luminosity function). 



4 EFFECTS OF STELLAR COLLISIONS ON COLORS OF GALACTIC CENTERS 

As mentioned in § 1, stellar collisions affect the luminosity properties and colors in galactic centers. Stellar collisions may 
disrupt one or both of the colliding stars to decrease the total luminosity, or the two colliding stars may merge to become 
a new star with luminosity properties or colors different from its parent stars. In this section, we will first give an general 
analysis of the effects of stellar collisions on color gradients, and then study these effects in nearby galactic centers, especially 
two galaxies in the Local Group: M32 and M31. 



4.1 Color gradients caused by stellar collisions 

Color gradients in a galaxy are usually caused by the variation of its stellar population with galactic radius. 

If there are no collisions in the stellar system, the surface brightness at color band C and at projected radius R, is given 

by: 



7o,c(-R,T age ) = J ' dz J dmda ^o(m, a,T agc )n(r,T age )Lc(m,a,T age ) (36) 
= Jo,c(R,T^ c )Rf (R). (37) 



Here r = \/R 2 + z 2 , £o{m, a, T agc ) is the expected stellar mass and radius function at the present time obtained by assuming 
an IMF and then only considering the effect of stellar evolution, Lc(m, (Z, T^ge) is the current luminosity emitted by a star 
with mass and radius (m, a) at color band C (in units of the solar C-band luminosity L©,c), ■/o,c(r,T , ag e) is the luminosity 
density at galactic radius r defined by 

Jo,c(r,T agc ) = /drnda ^o(m,o,T a ge)n(r,T ag e)Lc(m,o,Tag e ) (38) 

and fo is a factor reflecting the projection effect of the surface brightness defined by: 
f°° n(r,T a ~)dz 
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The color index of the surface brightness (at color bands Ci and C2) at projected radius R is given by: 
Mo,d(-R) - Ho,c 2 (R) = -2.5log[I 0}Cl (R)/Io,c 2 (R)] + (M , Cl - M Q ,c 2 ), 



(40) 



where M0,c ; is the solar absolute magnitude at color band d. We will call the surface brightness and the color in equations 
(36) and (40) the "original" surface brightness and color, that is, the brightness and color that would be present in the absence 
of stellar collisions. 

If two stars with masses and radii (rrii,ai) and (mj,aj) collide at time t' , we define Lij t c(t,t') (t ^ t') as the C-band 
luminosity of the collision product at time t. The luminosity and color of collision products depend on the types of colliding 
stars, their masses and radii, their relative velocity and impact parameter etc. If both of the two stars are disrupted by the 
collision, we have Lij t c(t, t') = 0; if the two stars merge to be a new star, Lij t c(t, t') is just the luminosity of the new merged 
star; and if one or both of the two stars survive the collision, L i j i c(t, t') represents the total luminosity of the surviving star(s) 
after the collision. Using the collision rate 1Zij(r,t') in equation (4), the change of the surface brightness at color band C 
caused by collisions of two stars with total mass larger than rriij is given by: 



-^coll,>mj 
C m t! 



J ir 



,o(-R,Tage) 

Ami 



r m tt i ff / ,Ta « 6 

/ dnij / ddi / ddj dz dt' lZij(r,t') 

Jmaxfmii-mj.O) J J J J () 



max(m^j — m t f ,0) J max(m^j — ,0) 

x [(Lij,c{Tnge, t' , r)) — Lc(?ni,ai,ragc) — Lc (rrij , aj , T a , 



(41) 



where {Lij t c(Ta, sc ,t' ,r)} is the average present luminosity of the collision product of two stars with masses and radii (mi,a,i) 
and (rrij,aj) colliding at time t' and at galactic radius r. We have the factor "1/2" in front of the integration in equation (41) 
for the same reason given for equation (8), and we have 
^ r m u r m tf r m tf r m a 

- / dirij / drrij — / dm, / drtij. (42) 

>/maxfmjj-m t f,0) ^ max(m ij — ,0) ^ max(m^j — m t f ,0) >< max(mij-mj,mj) 



The dependence of the average luminosity {Lij t c{T a&a ,t' ,r)} on r is due to the dependence of the collision outcome on the 
relative velocity of two colliding stars and the dependence of the relative velocity on r. In equation (41), we have ignored 
the radial motion of the collision products, because the radial motion of the collision products is expected to dilute the color 
gradients caused by collisions; thus, the color gradient obtained from equation (41) is an overestimate. We will mainly analyze 
the color gradient in the region with t co ii{r) ^> T age , where the stellar distribution and population [i.e., the stellar density 
n(r) and the stellar mass and radius function £(m, a)] evolve rather slowly, and we ignore collisions with collision products. 
Once the initial rapid stellar evolution is complete [massive stars have evolved into stellar remnants and the main-sequence 
stars have low masses (e.g., < IMq) with long lifetime in the main-sequence phase], in the slowly evolved region [where 
tcoii(r) » Tago], we have TZij(r,t) ~ TZij(r, 0) ~ TZ i:l (r, T agc ) in equation (41), i.e., 



-^coll,>mi. 
m tf 



:-(R,T at 



«/ ir 



drm 



max(mjj —m t f ,0) J max(m 

(iii,c(T'ago,t',r))dt' - Lc(mi,ai,T age )Tz 



drrij J dai J daj J dz IZij (r, T a( 



L 



Lc{m,j, Oj,T a| 



)Ta, 



Thus, the color index at projected radius R is given by 

., - 1 Io,c x (R) + -?coii,>o,Ci (R) 
MCi (R) - Hc 2 (R) = -2.5 log 



+ (M , Cl -M Q ,c 2 ), 



Io,C 2 (R) + Icoll,>0,C 2 (R) 

and the change of the color index caused by stellar collisions is given by: 
A/j co n,c 1 -c 2 (-R) = [/uci(-R) - vc 2 (R)] - [no,d{R) - im>,c 2 {R)]- 



(43) 



(44) 



(45) 



The contribution to the change of the color index from collisions between two stars with total mass rriij can be seen from the 
following analysis. According to equations (44) and (45), we have 



A/! co ii,Ci-c 2 {R) 



-2.5 log 



2.5 



/ 



1 + / d In nii 



-d/ co n, >mi . , Cl {R)/d\nmij 



Io, Cl (R) 



+ 2.5 log 



1 + 



In mi 



-d/ c< 



tj ,c 2 (R)/d In rriij 



In 10 



dlnm-i 



-d/coii, >m zj ,c\ (R)/d In rrii 
Io,cAR) 



h,c 2 (R) 
-dI co ii >mtj : c 2 (R)/d Intriij 



Io,c 2 (R) 



(46) 
(47) 
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dim 



■d7coii,>m 4j ,Ci(-R)/d In rriij 



Io,cAR) 



Afj, 



•coll , mi j ,C\ 



-c 2 (R), 



(48) 



where 



,>rriij ,Ci 

(R) 

din mi 



rn 



ij 



r m tt 



J max(m i j — m t f ,0) 



Ami I da; / da,j / dz IZij (r, T a 



/ {Lij,c{T& sc , t', r))dt' — Lc(m,i, a 4 , r ago )T a g C — Lc(mj,aj,T a , ge )T i 
Jo 

A/i co ll irrly ,c 1 -C 2 (R) = [A'coll, my,Ci(-R) — /icoll.m 
d (-R) — /Ucoll,my,C 2 (-R) 



Mcoll 



-2.5 log 



c 2 (-R)] - [Mo.Cj (-R) - Mo,c 2 (-R)], 
d/coii, (R) /d In my 



d/coii, > my , c 2 ( -R) /d In mi 



+ (M@, Cl - M ,c 2 ). 



(49) 
(50) 
(51) 



In equations (47) and (48), we use the relation ln(:r) — x — 1 if x — > 1 where Inx — A// co ii,my,Ci— c 2 /2.5. Equations 
(47) and (48) are valid for 7 co ii,>o,Ci {R)/Iq,c 1 (R) <S 1 and 7 co ii j>0 ,c 2 (-R)//o,c 2 (-R) <C 1; and for equation (48), we also 
need |A/i co n .my ,Ci — c 2 (7?) | 2.5. As seen from equation (48), for the collisions of two stars with the natural logarithm of 
their total mass in the range In — > In + d In , whether the change of colors caused by collisions is significant is 
determined by two factors: (i) the ratio of the change of surface brightness caused by collisions to the original surface brightness, 
dlnra» : j(d/ C oii,> mi:j ,Ci/dmmjj)/.Zo,C'; L , (ii) the difference between the color index of the change of surface brightness caused 
by collisions /icoii,m w ,Ci - Mcoii, m w ,c 2 and the original color index p ,Ci - fM),c 2 of the stellar system (i.e., A/u co ii 

In this paper, we focus on studying the effects on the color gradient in visible bands caused by collisions between main- 
sequence stars. We ignore the effects of collisions with post-main-sequence stars because the outcome (especially luminosity 
and color properties) of these collisions is not well understood and detailed studies of collisions with post-main-sequence stars 
(e.g., Bailey & Davies 1999) are beyond the scope of this work. We also ignore collisions of compact stellar remnants with 
main-sequence stars, post-main-sequence stars or even with other compact remnant stars, which usually form "exotic" objects 
and occur at low rates. We do not intend to study the color gradients in UV, f&r-UV or other invisible bands in realistic 
galaxies, because the luminosities in these bands are more sensitive to post-main-sequence stars (e.g., giants or horizontal 
branch stars) etc. 



4-1.1 Luminosity of collision products of main-sequence stars 

The outcome of collisions between main-sequence stars has been investigated by smoothed particle hydrodynamics simulations 
[e.g., see Sills et al. (1997, 2001) for collisions in globular clusters, and Freitag & Benz (2001) for collisions in galactic centers]. 
One of the possible outcomes is that the two colliding stars merge to form a new star, which is suggested to be one of the 
mechanisms to form blue stragglers [for other blue straggler formation mechanisms, see Leonard (1989), Stryker (1993) and 
Bailyn (1995) etc.]. Here, we simply call the merger products "blue stragglers" (BSs), whether or not they really correspond 
to the BSs found in observations. The BS is brighter and bluer than its parent stars for two reasons: it may have a larger 
mass, or the stellar envelope may have a higher helium content if its parent stars have evolved significantly before the collision 
and the central helium is mixed to the surface of the BS during the collision (e.g., Sills et al. 1997, 2001). The answers to 
the questions whether collisions of main-sequence stars can form BSs, how much mass is released from stars during stellar 
coalescence, and how much helium is mixed during the collisions, depend on the masses and radii of the colliding stars, their 
relative velocity and their impact parameter etc. Hence, these factors also determine what the luminosity properties of BSs 
formed by stellar collisions will be. Collisions of two stars with low relative velocity (e.g., < Vij in eq. 1) are more likely to lead 
to stellar coalescence and form BSs; and those with high relative velocity (e.g., at r < 10~ 3 pc in galaxies with M. ~ 10 6 M© 
or r < 1 pc in galaxies with M, ~ 10 9 M Q where the velocity dispersion is significantly higher than v esCj 0) are more likely to 
destroy stars, or the stars may pass through each other and both survive the collision (see Freitag & Benz 2001; Lai, Rasio 
& Shapiro 1993; Benz & Hills 1987, 1992). Below we will consider several extreme cases of collision products to study the 
effects of collisions on color gradients. None of the extreme cases are realistic, but they will simplify the complexity of the 
collision outcome and provide useful limits to the effects caused by stellar collisions. 

For the mass of the collision product, we consider two extreme cases: one is that collisions always destroy both of the 
colliding stars; the other is that collisions always produce BSs and the mass of the BS is the sum of the masses of its parent 
stars. 

We next consider the helium abundance of BSs. A star which has been on the main sequence for time r has created a 
mass of helium equal to r(L)/(ec 2 ), where (L) is the mean luminosity of the star during its main-sequence phase so far, and 
e = 0.007 is the efficiency of hydrogen fusion. We define Yq and Zq as the initial helium abundance and metallicity of the stars. 



© 0000 RAS, MNRAS 000, 000-000 



16 Q. Yu 



If a BS is formed by the collision of two main- sequence stars at time t and the chemical abundance is fully mixed (which is 
the extreme case of mixing), the resulting helium abundance Yij of the BS will be (see eq. 2 in Bailyn & Pinsonneault 1995): 

IWo+tro - 1 ^ {Li} ; {L >\ l-Yo), (52) 
rrii + rrij Lq 

where £i is the time in units of 10 10 yr, and and rrij are the masses of the two colliding stars. We will consider two 
extreme cases for the helium abundance of BSs: one is that BSs have the initial helium abundance of their parent stars 
Yb; and the other is that helium is completely mixed during collisions and BSs have the highest possible helium abundance 
Yij = O.ltio + Yo(l — O.liio) (eq. 52 for collisions of two main-sequence stars with solar mass). 

According to the above treatment, the luminosity of the collision product in equation (41) or (43) is bounded by the 
following three extreme cases: (i) Lij } c(t,t') = 0, (ii) 

Lij,c(t,t') = Lc(mi + m,j,t — t') (53) 

with helium abundance Y and metallicity Z , or (iii) Lij t c(t,t') = Lc(rrii + rrij,t — £') with helium abundance Y t j — 
O.ltio + Y (1 — O.ltio) and metallicity Z . Note that the merged BSs are assumed to be equilibrium main-sequence stars after 
collisions, because the BS should return to the main sequence on the short Kelvin- Helmholtz timescale (e.g., ~ 10 7 yr for stars 
with solar mass). In addition, a complete resetting of the nuclear reaction on the main sequence assumed in equation (53) 
(see the term "t — t'" ) is unlikely, especially for merged BSs with little helium mixing, which have a shorter main-sequence 
lifetime (Sills et al. 1997, 2001); and this assumption in equation (53) will give an upper limit to the effect of stellar collisions 
on color profiles. 



4.2 Color gradients in M32 

In this subsection, we use the dwarf elliptical galaxy M32 in the Local Group as an example to study the color gradients 
caused by stellar collisions in galactic centers. M32 (NGC 221) is useful for studying stellar collisions because (i) it is the 
closest elliptical galaxy, (ii) it is compact and its density is high, (iii) it has no visible features of dust or other anomalies 
(see details in § 4.2.1). M32 has the a shorter stellar collision timescale at HST resolution than any other nearby galaxy (see 
Fig. lb), and thus we expect that M32 will have the most easily detectable color gradients caused by stellar collisions (if any). 



4-2.1 Observations 

Both photometry (in the optical/infrared bands) and spectroscopy suggest that the center of M32 has a relatively homogeneous 
stellar population without strong gradients in age or metallicity. No evidence of features such as an inner disk, dust, or any 
other structure is visible in the HST WFPC2 optical images and the Near Infrared Camera and Multi-Object Spectrometer 
(NICMOS) images (Lauer et al. 1998; Corbin, O'Neil & Rieke 2001). A Nuker-law fit (eq. 19) to the V-band surface 
brightness profile at r < 1" gives a = 1.39±0.82, = 1.47 ±0.16, 7 = 0.46 ±0.14, I b corresponding to ^ b = 12.91±0.31, and 
r b = 0.47" ± 0.15" (Lauer et al. 1998). The central cusp of the surface brightness profile of M32 continues to rise into the 
HST resolution limit; thus there is only a lower limit on the central density po > 10 7 Mq pc~ 3 (Lauer et al. 1998). The BH 
mass of M32 is about (2.5 ± 0.5) x 10 Mq and the stellar /-band mass-to-light ratio (determined from dynamical models) is 
M/Li = (1.85±0.15) M /L o ,j (Verolme et al. 2002). The fit to the M32 color profiles for 0.1" < r < 10" gives an essentially 
flat color index as a function of projected radius R (Lauer et al. 1998): 

Mobs.v - Atobs,/ = (1-237 ± 0.002) - (0.009 ± 0.005) log(i?/l"). (54) 
and 

Hobs,u - Mobs.v = (1.216 ± 0.004) - (0.023 ± 0.008) log(i?/l"), (55) 

where /i D b s is the observational surface brightness in magarcsec~ 2 . The absorption line indexes in the optical bands also show 
no radial gradients from 1" to ~ 10" (del Burgo et al. 2001). HST NICMOS imaging of the core of M32 shows that the 
near-infrared surface brightness profiles can be fitted by the same form as that fitted to the optical images (Corbin, O'Neil & 
Rieke 2001) and no strong gradients are found in the infrared color profiles (Peletier 1993; Corbin, O'Neil & Rieke 2001). 
The spatial distribution of the brightest stars in the infrared also follows the optical light profile outside 2" of the nucleus of 
M32 (Davidge et al. 2000). The absence of gradients in both colors and absorption line indexes indicates that the stars are 
coeval or that the stars have been mixed considerably in the inner region of M32. By fitting the mean values of the spectral 
indexes and colors for the inner regions of M32 to the stellar population synthesis model in Vazdekis et al. (1996), del Burgo 
et al. (2001) find that the population in the core of M32 can be modeled as a single stellar population of an intermediate age 
(~ 4Gyr) and solar abundance. Note that the IMFs used in the model of Vazdekis et al. (1996) are not exactly the same as 
equations (26) and (27) (e.g., they have different lower and upper limits of the stellar mass range etc.). In this paper, we will 
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focus on the results obtained by using the Salpeter IMFs in equation (26), but our conclusions will not be significantly affected 
if we use the IMFs in Vazdekis et al. (1996) or equation (27). We will study whether collisions between main-sequence stars 
can produce a color gradient in M32 in visible bands. 



4-2.2 Expected color gradients caused by collisions between main-sequence stars 

Assume that the center of M32 is composed of a single stellar population of an intermediate age (~ 4 Gyr) with a Salpeter 
IMF (eq. 26) and solar abundance (see § 4.2.1), and that the stellar population evolves as described in § 3.4. In such a system, 
the turn-off stars have a mass m t f ~ 1.3 M© (set £ms = 4 Gyr in eq. 28). 

In § 4.1.1, we have described three extreme cases of the luminosity of collision products Lij.c(t l ', t). Using equations (36), 
(40)-(45), and the luminosity of main-sequence stars Lc(m,t) obtained from the Padova stellar evolutionary tracks (Girardi 
et al. 2000) , we do the numerical calculation for cases (i) and (ii) and find that stellar collisions cannot cause observable color 
gradients in color indexes fj,u — fJ-v and fj,y — f^i in the center of M32 (i.e., A/u co ii,(/-v < 0.02 mag and A/v>n,v-i < 0.02 mag 
at R > 0.1", see Fig. 10 below). In case (iii), the evolutionary tracks of stars with the required chemical abundance is not 
provided in Girardi et al. (2000) or any other published literature that we are aware of. However, we shall argue that the 
difference in the chemical abundance in case (iii) will not significantly affect our conclusions obtained in cases (ii). 

The above numerical results (i.e., no observable color gradients are expected to be caused by collisions between main- 
sequence stars in M32) may be understood in the following semi-analytic way. We consider the case in which stellar collisions 
always lead to stellar coalescence, i.e., case (ii). Using equation (53), we have J ago Lij t c(Ta, ge , t') At' = J as ° Lc(rrii + 
n r ij,T a .g e — t') dt' = J Tagc Lc(rrii + m,j,t) dt, and the surface brightness contributed by the BSs which are formed by collisions 
of two stars with total mass higher than rriij is given by (cf., eq. 43): 

Icoll,>m i:j ,c(R, Tago) 

= - / dm; / drrij / da* / do,- / dz TZij (r, T ago ) 

Jmax(m i j-m t f,0) J maxfrn^j — ,0) J J J—oo 

x / Lc(rrn +m,j,t) dt (56) 



— J r ra?i,>m i:Ji c(-R,T' a g C )i?/coii(-R) (57) 

where J^fi^my ,c(r, T^ gc ) is the luminosity density contributed by the BSs which are formed by collisions of two stars with 
total mass higher than rriij at galactic radius r: 

Jco\\,>ma ,c(r : Tage) 



= - drrii / drrij / dat / da., 7^(r,T a( 

Jmax(m,i-m t f,(l) </ maxfm,- t — ,0) J J 

"/ 

JO 



(m^j— m t f,0) J maxjmj. 

T'agc 

Lc{nii + nij,t') dt' . (58) 



Here / co ii(-R) is a factor reflecting the projection effect of the surface brightness. According to equations (l)-(4) and (56)-(58), 
we approximately have 

J-°L n 2 (r,r agc )a(r)[l + v 2 asc , G /6a 2 (r)] dz 
/coll(i?) ^ n*(R,T^)o-(R)[l + vi Ct J(>o-HR)]R ' (59) 

where v t j ~ tw,o, (^rci) a(r) and (v re i) (v^ 1 ) 1 — (« 2 i) = 6<7 2 (r) [/ co ii(-R) is not sensitive to the factor "6" in front of 
cr 2 (r)] are assumed (see eq. 1). Thus, as in the derivation in equations (44)-(51), we may obtain the change of the color index 
caused by BSs given by: 

A bs , m /coii(fl) f -dJ^ Cl (i?)/dlnm i3 

AMcoii,Ci-c 2 (-R) ^ T7d\ / dlnm H T — 7m A^ coll , m Cl -c 2 (i?) (60) 

Jo{R) J Jo,cA H ) 

where 

dJ^ll^mij ,Ci (R) m, 



Jo 



,t r r 

drrii / dai / daj TZij(r,T age ) 

max(mjj -m t f,0) J J 

L c {m tJ ,t) dt, (61) 



A/icoU.my.d-CaCR) = [Vcoll, m^.diR) - AWl, mi j ,C 2 {R)\ ~ [A»0,Ci (-R) ~ H0,C 2 (R)], (62) 
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McoU.my ,Ci (R) ~ l^oll,m, ij ,C 2 {R) 

r dJ c B f 1 , >mi3 ,c 1 (fi)/dlnm I3 



-2.5 log 
-2.5 log 



L Cl (m tJ ,t) dt 
Lc 2 (m ijy t) dt 



+ (M , Cl -M ,c 2 ) 
+ (M Q , Cl -M ,c 2 ). 



(63) 



Equation (60) is valid for [/ co ii(fl) JcoTi^o.Ci (R)]/[fo(R) J , Cl (R)] < 1, [/coii(-R) ^ co Ti,>o,c 2 (-R)]/ [/o(-R) Jo,c 2 (-R)] < 1 and 
l^/ u «>n."iy.c r i-c 2 (^)| "C 2.5. Equation (60) shows that for the collisions of two stars with total mass in the range In rriij — ► 
In iriij + dlnm^-, whether the effect of BSs on the color gradient is significant at projected radius R is controlled by three 
factors: (i) the ratio of the luminosity density d In rriij (dj^ } > mij ,Ci /dlnrriij) / Jo,Ci at galactic radius r = R; (ii) the difference 
between the color of BSs and the original color of galactic centers, A/i^ 11}rnij} c 1 -c 2 ', and (iii) the projection effect / C oii//o 
(see eqs. 39 and 59). 

Below, we will give a quantitative estimate to the three factors in equation (60). If the center of M32 is composed of 
identical stars with solar mass and radius, we may use the surface brightness profile of M32 and its mass-to-light ratio T 
determined from dynamical models to get its stellar mass density or stellar number density n (r). In such a system, we 
define the luminosity density contributed by BSs formed by collisions between solar-type stars as Jf^\ t Q,c(r); and according 
to equation (43), we have 



J. 



coll,0 



,c(r) 



[n Q (r)/2£ coll , (r)]/ o ago Lc(2M ) dt T c f MS (lM ) 



Jo,c(r) 



0.2 



n (r)L c (M )/T c 

10" yr tcoii,o(r = 0-l") T c 
t co ii, (r = 0.1") tcoii, (r) 2 



icoii, (r) 



(64) 



where J Tag ° Lc(2Mq) dt ~ Lc(l M )iMs(l M ) is assumed (i.e. the time-integrated luminosity of main-sequence stars is 
assumed to have / () T ° SC Lc(m,t) dt oc m for m > m t f), and £ms(1M ) = 10 10 yr. If the stars have a distribution of masses 
and radii, according to equation (58), we have the ratio of the luminosity density contributed by BSs as follows (similarly as 
eqs. 31, 33 and 35): 



tBS 

- / coll,>m i 



,c(r) 



rBS 

" / coll, 
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[(a l + a,)/(2R )] 2 

[(a i +a 3 -)/(2R )][(m i + m j )/(2M )] 



drrij H(mi, 0)H(mj, 0) 



J o Tago L c {mi+m jy t) dt 
J Tage L c (2M e ,t) dt 



if o(r) > iw, 
if cr(r) < w esc , 



We assume the following relation for the luminosity of main-sequence stars: 



/-"age . 

I Lc(m, t) dt oc < 



m if m ^ rriti 
m VL if m < mtf 



(65) 



(66) 



(i.e., the luminosity of main-sequence stars is assumed to be Lc(m,t) oc m VL ; t)l = 3.5). Using equations (65) and (66), we 
show in Figure 7 the ratio of the luminosity density J^\ y>mjj ,c / Jfoii,Q,c (solid lines) and its derivative \djf^ ll}>m .. ,c /dlnm^l 
/Jc1i,q,c (dotted lines) as a function of stellar mass rriij. We show the results for both of the two cases cr(r) ^ iw.o and 
cr(r) <C «csc, in equation (65), and we have 



Jcol,>o,cjr) 



0.3. 



(67) 



The peak of the derivative IdJ^ t>mij ,c /din niij\/ Qi c is located at rriij ~ m t c, which represents that most of the 
luminosity contributed by BSs comes from collisions between two main-sequence stars with total mass rriij ~ m t f. 

To obtain the difference between the color of BSs and the original colors of the center in M32 in the absence of stellar 
collisions (i.e., A/i^ ljmi:ji Ci-c 2 in eq. 62), we use the Padova stellar evolutionary track (Girardi et al. 2000) to obtain the 
color indexes of time-integrated luminosity J^ s ° Lc (rriij, t) dt of stars during their main-sequence phases as a function of 
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stellar mass my, Ad (my) -Ac 2 (mij), defined by (cf., eq. 63) 

L Cl (m l3 ,t) dt 



Ad (my) - Ac 2 (mij) = -2.5 log 



Lc 2 (m t „t) dt 



+ (M , Cl -M ,c 2 ) (68) 



(in the integration / Tas ° Lc(mij,i) dt, even if we include the luminosity during the post-main-sequence phase, our conclusion 
will not be affected). Combining equations (63) and (68), we have 

Mcon.my ,Ci - McoU,m w ,c 2 = Ad {rriij) - Ac 2 (mij). (69) 

The A(/(m»j) — Ay(my) and Ay (my) — A/ (my) for stars with solar abundance are shown as a function of stellar mass rriij in 
Figure 8. The horizontal lines in Figure 8 give the observed color indexes of fj, bs,u — Mobs , v (thin solid line) and piobs , v — Mobs ,/ 
(thin dashed line) in the center of M32 at 7? = 1" (eqs. 54 and 55). As seen from Figure 8, for color index Au (rriij) — Ay (my) 
(thick solid line), the stars with mass m tJ > 0.9 Mq are bluer than the center of M32; and for color index Ay (my) — Ai (my) 
(thick dashed line), the stars with mass my > 0.7 M Q are bluer than the galactic center of M32. In Figure 8, at my ~ m t f, 
where the peak of |d7^,f li>m ,c/dln my |// C oii,o,c is located, we have 

A//coU,my,i/-v = [A(/(my) - Ay (my)] - [^obs,[7 - Atobs.y] - -0.7 (70) 
and 

A/icoU.my.v-/ = [Ay (my) - A/ (my)] - [^ bs,y - Mobs,/] ^ -0.6. (71) 

Using equations (39) and (59), we obtain the projection factors fo{R) and f co ii(R) of M32 and their ratio f co ii(R) / fo(R) as 
a function of projected radius R, which are shown in Figure 9. The results in the region not resolved by the HST (2R < 0.1") are 
obtained by extrapolating the Nuker-law fitted surface brightness profile inwards. Figure 9 shows that f co n(R = 0.1")//o(7? = 
0.1") weakly depends on 7? and 

/ coll (7? = 0.1") _ 

MR = o.i") 

Using equation (60), we have the following relation for the change of the color index caused by stellar collisions: 

I A BS >mi . /coll (7?) Joon>0,Ci (-R) | A BS | /- Q -, 

lAMcoii.Ci-cM-R)! < J^ftj J 0C (R) l A ^oii,m tf ,Ci-c 2 |, (73) 

where we use the value of A/i c 3 n,m i:j ,Ci-C2 , at my ~ m t f since the peak of |d/^,fi i>m4 ./dlnmy|/7o is located at my ~ m t f 
(note that ^7^ >mi ./dlnmy|/7o decrease sharply from my ~ m t c to both high-mass and low-mass ends in Figure 7 and 
Ad (rriij) — Ac 2 (my) change mildly in the whole mass range in Figure 8). Using equations (64), (67), and (70)-(72) in equation 
(73), we have 



A^iuu-v(R) 



°> i a bs fri\ >> -0-03 mag 



tcoll,0(7£) 



tcdU,©(fl = 0.l") 



(74) 



where the collision timescale t co ii,0 (7? = 0.1") ~ 10 11 yr (see Fig. 1) for M32 is used. Inequality (74) gives an upper limit of 
the blueing caused by collisions between main-sequence stars, which shows that it is almost impossible for stellar collisions to 
cause observable blueing at R > 0.1" in the center of M32. But inequality (74) does not exclude the possibility that stellar 
collisions can cause a color difference > 0.02 mag in the much inner region at 7? < 0.1". Note that inequality (74) is invalid in 
the region [e.g., 7? < 0.01" with t C oii(7?) < T agc ] where the approximation in equations (48) and (60) is invalid. 

In the assumed extreme case (iii) on the luminosity of collision products, BSs have a different Helium abundance from 
case (ii). But the difference in Helium abundance will not significantly affect the luminosity and colors of BSs or the terms 
jf<A\,>m i:j ,c(R)/Jo,c(R) and A^aii,m ij ,c 1 -c 2 (R) (cf-, eq. 60), and hence will not affect our conclusions for the following 
reasons. In cases (ii) and (iii), the mean molecular weights /Kw = (2X + 0.75F + 0.5^) _1 are about 0.62 and 0.63 (Helium 
abundance Y = 0.28 and 0.31), respectively. For main-sequence stars with same mass, the increase in the mean molecular 
weights can only increase the luminosity by about eight percent (L oc fx%, see eq. 20.20 in Kippenhahn & Weigert 1990). The 
radius of the stars (a* oc Mw' 6 , see eq. 20.21 in Kippenhahn & Weigert 1990) depends weakly on the mean molecular weight. 
Thus, the increase in the mean molecular weight and stellar radius will cause the increase of the effective temperature of the 
star (T e ff oc L/R 2 oc Mw 8 ) by about five percent, which may make main-sequence stars with mass around m t { (with effective 
temperatures T e g ~ 6x 10 3 K) bluer only by less than 0.1 mag in Au — Ay and Ay — A/ (or A/ttfon, m tt ,u-v and AMc n,m tf ,v-i)- 
Therefore, the results obtained for case (iii) (considering the change of the luminosity and colors of BSs caused by helium 
mixing) should not be significantly different from those obtained for case (ii). 

Inequality (74) only provides the effect on color indexes caused by BSs. Including the effect of removing the luminosity of 
parent stars from the galactic center [i.e., "the term —Lc(m i ,ai,T l> ,g C ) — Lc(m,j,aj,T agc )" in eq. 41] will not make the galactic 
center become bluer because the decrease of the luminosity caused by destruction of main-sequence stars comes mainly from 
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destruction of high-mass stars (just as most of the mass density involved in collisions come from destruction of high-mass 
stars shown in Fig. 7), which are usually bluer than the center of M32 (cf., Fig. 8). If collisions of main-sequence stars are 
assumed to always lead to destruction of main-sequence stars and no BSs are formed [i.e. case (i)], using the above similar 
analysis for the effects by BSs, we find that there will be no observable reddening at R > 0.1" of M32, either. 

Figure 10 shows our numerical results for the color indexes fiv{R) — l-ii(R) and Hu{R) — ilv(R) in M32, obtained by 
using equations (36), (40)-(45) and considering the luminosity of collision products in case (i) and case (ii) (dashed line and 
solid line) as described in § 4.1.1. We use the observed color indexes (eqs. 54 and 55) as the original color index /xo.Ci — Mo,c 2 • 
Figure 10 shows that collisions of main-sequence stars cannot cause observable color changes in the visible bands at R > 0.1" 
in M32, which is consistent with the observation that no color gradients are seen in M32 by HST; and at even smaller radii, 
stellar collisions are likely to cause the color difference larger than > 0.02 mag (e.g., at most 0.08 mag at R — 0.02"). The 
above analysis and numerical results are only valid if t co n(R) > T age , In the region with £ C oii(-R) J; Tage (e.g., R < 0.01"), it 
is possible that most of the main-sequence stars have experienced collisions. If we extrapolate the observed V-band surface 
brightness profile inward and assume that the surface brightness in the region with £ C oii,0(-R) < T ase has the color of BSs 
described in equations (70) and (71), and then use the above numerical results of the surface brightness in the region with 
i C oii,©(r) > Tagc, we may obtain the averaged color index within a given galactic projected radius R. We find for case (ii) that 
the averaged blueing within R < 0.1" due to collisions is not larger than 0.06 mag in fiu — fiy and 0.02 mag in /iv — /i/, and 
the averaged blueing within R < 0.05" due to collisions is not larger than 0.16 mag in fiu — l-iv and 0.06 mag in fiv — la- 
in Lauer et al. (1998), stellar collisions are claimed to be capable of causing observable blueing in the center of M32. 
However, Lauer et al. did not account carefully for distribution of stellar masses and radii, projection effects, and the difference 
between the color of BSs formed by stellar collisions and the original color of the center of M32. 

Based on the above results for M32, we expect that stellar collisions cannot cause any observable color gradients at 
R > 0.1" in visible bands in the centers of M31 and other nearby galaxies. This result may be understood as follows: 
these systems have longer collision timescales than M32 (cf., Figure, lb or Table 6 in Lauer et al. 1998); compared to 
the values obtained for M32, the factor |A/icoii,m t f ,Ci-c 2 1 (the difference between the color of BSs with mass m t f and the 
original color of galactic centers) in equation (73) should not change significantly; and the ratios of the projection factors 
fcoii(R <C r^)/fo(R <S fb) for core galaxies with small 7 may become even smaller, and those for power-law galaxies do not 
significantly increase. Multicolor HST WFPC2 images show that the center of M31 appears bluer in a 0.14" x 0.14" box than 
the rest of the nucleus and the surrounding bulge (1" < R < 5"), by ~ 0.1 mag in fiv — l^i (see Table 2 in Lauer et al. 1998). 
Our estimate shows that the averaged blueing is not larger than 0.03 mag in fiv — Hi within R < 0.14"/2 even for M32, and 
hence the ~ 0.1 mag blueing in M31 is unlikely to be caused by stellar collisions. It is therefore worthwhile to investigate other 
mechanisms which may cause the blueing in the center of M31. 



5 CONCLUSIONS 

We have studied the effects of stellar collisions in realistic galactic centers, focusing on feeding massive BHs and color gradients. 

We find that the mass involved in stellar collisions is not sufficient to provide the present BH mass in realistic galactic 
centers, especially in bright core galaxies. Similarly, Magorrian & Tremaine (1999) showed that stellar tidal disruption rates 
cannot contribute significantly to the present BH mass in galaxies brighter than ~ 10 9 L Q . These two results together discredit 
two proposed mechanisms (stellar collisions and tidal disruption of stars by central BHs, e.g., Frank 1978) for the growth 
of massive BHs in the centers in bright galaxies. The raw material for BH growth must therefore come from other sources, 
for example, the mass released by stellar evolution in the initial ~ 1 Gyr of the galaxy's lifetime or gas that sinks to the 
galactic centers in a merger (see simulations on galaxy encounters in Barnes & Hernquist 1996). For both of the above 
mechanisms, massive BHs may form and grow rapidly, early in the history of the universe (e.g., at redshift z ~2-3, where 
galaxy encounters mostly occur in the hierarchical galaxy formation model). This scenario is consistent with the argument 
made by Yu & Tremaine (2002) that growth of high-mass BHs (> 10 8 M Q ) comes mainly from accretion during optically 
bright QSO phases (which is obtained by studying the relation between the local BH mass function and the QSO luminosity 
function) . 

We have analyzed how the color of a stellar system is affected by collisions of main-sequence stars. We account for the 
effects of collisions of main-sequence stars on color indexes at projected radius R due to three factors (see eq. 60): (i) the ratio 
of the change of luminosity density caused by stellar collisions to the original luminosity density of galactic centers at galactic 
radius r = R; (ii) the difference between the colors of BSs and destroyed stars and the original color of galactic centers; and 
(iii) the projection effect between the surface brightness and the luminosity density. We find that collisions between main- 
sequence stars cannot cause observable color gradients in the visible bands at R > 0.1" in M32, M31 and other nearby galactic 
centers. This result is consistent with the lack of an observed color gradient in M32 at HST resolution. At even smaller radii, 
collisions between main-sequence stars are likely to cause the color difference larger than > 0.02 mag (e.g., at most 0.08 mag 
at R = 0.02"). The averaged blueing due to stellar collisions in the region R < 0.1" of M32 should not be larger than 0.06 mag 
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in nu — fiv and 0.02 mag in fj,y — Hi- The observed blueing in the center of the galaxy M31 (in a 0.14" x 0.14" box; Lauer et 
al. 1998) must be caused by some mechanisms other than collisions between main-sequence stars. 

This project was suggested by Scott Tremaine, my thesis advisor; and I am deeply indebted to him for his advice. I also 
thank Marc Freitag, Jeremy Goodman, Alison Sills and David Spergel for helpful discussions. Support for Proposal number 
HST-AR-09513.01-A was provided by NASA through a grant from the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy, Incorporated, under NASA contract NAS5-26555. Support was 
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Figure 4. The ratio of collision timcscales t co n i/t co \\ q as a function of stellar mass m;. The quantity t co u j (sec cq. 5 or 31) is the 
collision timcscale of stars with mass rrii obtained by assuming that the stellar system has a distribution of stellar masses and radii 
and ignoring collisions with post-main-sequence stars and stellar remnants; and t co \\ q (see Fig. 1) is the collision timescale of stars by 
assuming that stars in the stellar system have solar mass and radius. Note that t co ii j at m; = Mq is not necessarily equal to t co n q. The 
dotted curves represent the results for <x(r) <C tw,0, and the solid curves represent the results for cr(r) ^ u C sc,©- The curves ending at 
0.08 are obtained by assuming the Salpctcr IMF (eq. 26), and those ending at 0.01 Mq are obtained by assuming the multi-power-law 
IMF (eq. 27). The turn-off stars are assumed to have a mass m t f = 1 Mq. 
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Figure 5. Panel(a): The ratio of P^ifj/(2/5™^g) as a function of stellar mass m; (cf., eqs. 32 and 33). Panel (b): The ratio of 
P'ccM? / f'ccM Q as a f unc tion of stellar mass m; (cf., eqs. 34 and 35). The quantity p™* x 4 dlnra; is the mass density of main-sequence stars 
involved in collisions with main-sequence stars with mass in the range lnm, -> lnm; + dlnm^ per unit time; the quantity p' c ^fd\nmi 
is the mass density of main-sequence stars with mass in the range In m; —* In m.; + d In mi involved in collisions per unit time; and the 
quantity P™^ X q is the mass density of main-sequence stars involved in collisions per unit time that is obtained by assuming that all the 
stars in the stellar system have solar mass and radius. The ratios p™^ x 4 / '(2p™^ q) and P^JjiT Ip'Saiq are no * necessarn y 1 a * mi = Mq. 
The curves ending at 0.08 Mq are obtained by assuming the Salpeter IMF (eq. 26) and those ending at 0.01 Mq are obtained by assuming 
the multi-power-law IMF (eq. 27). The dotted curves represent the results for cr(r) <C v CBCi q, and the solid curves represent the results 
for a(r) > u csc ,©. 
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Figure 6. Stellar mass tidally disrupted or swallowed whole by central BHs M t id,Q(tHubble) versus stellar mass involved in collisions 
M™j5 x Q (tHubble) °ver a Hubble time (see Fig. 3a). All stars are assumed to have the solar mass and radius. Dotted lines are the reference 
lines for JWy^/M™^ = 0.1, 1, 10. Solid circles represent core galaxies and open circles represent power-law galaxies. My^Q and A^™J\ X q 
are generally comparable. See discussion in § 3.5. 
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Figure 7. The ratio of the luminosity density contributed by BSs formed by collisions between main-sequence stars, 
■^«>11 >m -.c/^Sfl 0.C' as a f unc tion of stellar mass rriij, in M32 (solid lines for a(r) <C v CSCi q or a(r) S> f C sc,©, cf., eq. 65). The 
■^coll >m >C re P resen ts the luminosity density contributed by BSs which are formed by collisions of two main-scqucncc stars with total 

mass higher than mass rriij, and M32 is assumed to have a Salpeter IMF (eq. 26) and age T ago = 4 X 10 9 yr (see § 4.2.1). The t/^jjj g,c 
represents the luminosity density contributed by BSs formed by collisions by assuming that the center of M32 is composed of identical 
stars with solar mass and radius. The dotted lines represent |dJ^ >m ,c/d lnmjj |/ J^fj g,c; an d their peaks are located at rriij — m ti, 
which shows that most of the surface brightness contributed by BSs comes from collisions between two stars with total mass around m t {- 
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Figure 8. The color indexes Ajj — Av (thick solid line) and Ay — ^1 (thick dashed line) of the time-integrated luminosity (Jq^ s " Lc dt) 
of stars during their main-sequence phases, as a function of stellar mass mij (eq. 68). The results are based on the stellar evolutionary 
tracks of Girardi ct al. (2000). The stars have solar abundance. The horizontal lines give the observed color indexes, l^ b a ,U ~ Mobs.v 
(thin solid line) and /Li bs,V — Mobs,/ (thin dashed line) of the surface brightness of M32 at R = 1" (eqs. 54 and 55). The vertical dotted 
line represents stellar mass m\j = 1.3 Mq, which is the turn-off mass at age T age ~ 4Gyr (see eq. 28). 




Figure 9. The projection factors fo(R) (long dashed line), / co ii(-R) (short dashed line) and their ratios fcoll(R)/fo(R) (solid line) as a 
function of projected radius R in M32 (see eqs. 39 and 59). The vertical dotted line represents R = 0.1". The results in the region not 
resolved by the HST (2R < 0.1") are obtained by extrapolating the Nuker-law-fitted surface brightness profile inwards. 
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Figure 10. The estimated color indexes fiy — fii and /ijj — fiy (eq. 43) in the center of M32, as a function of projected radius R. 
Equations (36) and (40)-(45) are used in the calculation. The dashed lines represent the color gradients obtained by assuming that 
collisions always lead to destruction of both colliding stars and that no BSs are formed, [case (i) of § 4.1.1]; and the solid lines show the 
results obtained for case (ii), in which the colliding stars are assumed to form a BS with solar abundance. The left vertical dot-dashed 
line marks the HST resolution (R = 0.0455") and the right vertical dot-dashed line marks R = 0.1". The solid squares are the color 
indexes obtained from observations (see Fig. 18 in Lauer et al. 1998); and the dotted lines represent a fit to the observed color index in 
the region 0.1" < R < 10" (eqs. 54 and 55), which is used as the initial state in our calculations (i.e., /jq,Ci — M0,C 2 m e 1- 40 and 45). 
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